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Abstract 

We show that all scaling quantum graphs are explicitly integrable, i.e. any 

one of their spectral eigenvalues E n is computable analytically, explicitly, 

and individually for any given n. This is surprising, since quantum graphs 

are excellent models of quantum chaos [see, e.g., T. Kottos and H. Schanz, 

Physica E 9, 523 (2001)]. 
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Graphs are networks of bonds and vertices. Figure 1 shows two examples: a three- 
bond four- vertex star graph (Fig. la) and a three-bond four- vertex linear graph (Fig. lb). 
A quantum particle moving on the graph turns the graph into a quantum graph [1]. If 
the quantum particle moves freely on the graph, subjected only to flux conservation at its 
vertices, we call it a standard quantum graph. This is the type of quantum graph most 
frequently studied in the literature [1-5]. A larger class of quantum graphs, including the 
standard quantum graphs, are dressed quantum graphs [6]. A dressed quantum graph has 
potentials on its bonds and 5-functions on its vertices. The potentials on its bonds are 
essentially arbitrary as long as they do not introduce turning points on the bonds. But even 
this case can be dealt with trivially by re-defining the topology of the graph. 

An important subset of dressed quantum graphs are scaling quantum graphs [7-9]. In 
this case the graph bonds are dressed with scaling potentials and the graph vertices are 
dressed with scaling ^-functions. A scaling potential is one whose strength Vo scales with 
the energy E of the quantum particle according to V = XE, where A is a constant. The 
strength of a scaling 5-function scales with k = y/E. Scaling potentials and 5-functions are 
a natural choice to consider. On the one hand they frequently occur in physical systems 
[10-15], on the other hand they are mathematically convenient, since they allow studying a 
quantum system without causing phase-space metamorphoses [16] in the underlying classical 
system. It has been pointed out before [17] that this is the most natural way of studying 
quantum systems, in particular quantum chaos [18,19]. Since quantum graphs are popular 
and successful models of quantum chaos [1-5], it may come as a surprise that the energy 
spectrum E n , n = 1, 2, ... of all scaling quantum graphs is explicitly and analytically solvable 
in the form E n — . . ., involving only known quantities on the right-hand side. In many cases 
the solutions can be stated in closed analytical form. 

The spectral function g^\k) of a general scaling, dressed quantum graph is of the form 

[7] 

g (0 \k) = cos(S k - 7r 7o ) - ]T af cos{Sjk - vr 7j ), (1) 
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where £0 > is the total reduced action length of the graph [7,8], < Sj < S are certain 
combinations of the reduced bond actions [7,8], iV is the number of action combinations in 
(1), 70, 7j are constant phases and af^ are constant amplitudes. The spectrum E n of the 
quantum graph is obtained by solving the spectral equation 

/WHO, n = l,2,... (2) 

via E n = (k^) 2 . For the purposes of this paper we are only interested in the positive 
solutions of (2), and obtain a well-defined counting index n by defining to be the first 
positive root of (2). As a first step toward the solution of the general problem, it was shown 
in [7-9] that (2) can be solved explicitly in the form k^ = ... if the regularity condition 

E i4 0) i < 1 (3) 

3=1 

is fulfilled. In order to substantiate our claim that (2) is solvable explicitly for all scaling 
quantum graphs, we have to show that (2) is solvable explicitly even if (3) is not fulfilled. 

Before we turn our attention to the general case, we introduce our methods with the help 
of a simple example. Let us consider a scaling quantum graph derived from the three-bond 
star graph shown in Fig. la by putting the scaling potentials Vi(E) = XiE, < A; < 1 on 
its three bonds of length Li, I = 1,2,3, require the "Kirchhoff-type" [5] flux conservation 
condition X)f=i dipi/dxi = at its central vertex (ipi is the quantum wave function on bond 
number / of the graph and X\ is the coordinate on bond number /) and require Dirichlet 
boundary conditions on its three dead-end vertices. The spectral equation is of the form (1) 
with iV = 3, 71 = 72 = 73 = and 
So = Oil + a 2 + a 3 , Si = — Oil + a 2 + a 3 , S 2 — a± — a 2 + a 3 , S 3 = <y,\ + a 2 — a 3 , 

(4) 

(o) _ -A +P2 + P3 (o) _ P1-P2 + fh (o) _ P1+P2- fh , s 

where 
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= ftL h A = y/l - A <> i = 1,2,3. (6) 



The amplitudes in (5) do not fulfill the regularity condition (3). In some cases Y%=i \ a f^\ = 1 
(for instance for af^ > 0, j = 1,2,3), and in many cases Y^=i \ a f^\ > 1, which strongly 
violates the regularity condition (3). Since the methods and techniques presented in [7-9] for 
obtaining the spectrum of a graph explicitly depend crucially on (3), it seems that completely 
different methods have to be developed for general graphs, such as the three-bond star graph 
of Fig. la, which do not fulfill (3). There is, however, a way to reduce (1) to a form that 
allows to bring the powerful theory of regular quantum graphs [7-9] to bear. In order to 
motivate and to illustrate this method, let us study the case a\ = 1, a 2 = 7, 0:3 = 11, 
p 1 = 1/10, (3 2 = 1/5, p 3 = 1/2. In this case = 3/4, a { 2 0) = 1/2, 4° } = -1/4, S = 19, 
Si — 17, S2 — 5, S3 = —3 and the spectral equation is given by 

g (°\k) = cos(19A;) - ^ cos(17A;) - - cos(5A;) + - cos(3A;). (7) 

4r Zi t: 

Since + (a^l + l ^! = 3/2 > 1, this quantum graph is certainly not regular. But let 
us look at the first derivative of (7). Dividing by So, we obtain 

3 

^(jfe) = cos[S k + tt/2] - ]T af co^Sjk + tt/2] = 

3=1 

515 3 
- sin(19fc) + — sia(17fc) + — sin(5A;) - — sin(3A;). (8) 
76 38 76 

This time we have Y%=i \ a f^\ — 16/19 < 1 and therefore, since (8) is precisely of the form 
(1) and satisfies (3), it can be solved explicitly using the methods of [7-9]. In particular 
it was shown in [7-9] that root number n of a spectral equation that satisfies (3), such as 
(8), is found in the root interval [k n -i,k n ], where k n are the root separators [7-9]. It was 
also shown in [7-9] that the location of the root separators is entirely controlled by the local 
extrema of the trigonometric function with the largest action argument. Thus, in our case, 
the root separators of (8) are given by k n = (2n + l)7r/38. Since according to [7-9] root 
number n and only root number n is located in the interval [k n -i, k n ], we can now compute 
all roots of (8) explicitly and individually according to 



1.(1) 



&n— 1 



dgU(k) 



dk 



5{g {1 \k))dk. 



(9) 



In [20] we show that because of the hermiticity of the spectral eigenvalue problem on quantum 
graphs the locations of the local extrema of g <y0 \k) are separators for the roots of g^\k). The 
location of the local extrema of g ( -°\k), however, are given by the zeros of g'" l \k), which, up 
to constants, is the derivative of g^(k). Therefore, using the roots k£\ explicitly computed 
in (9), as the root separators of (7), we obtain, again explicitly and individually, 



fc(°) - f 



k 



dg<®(k) 



dk 



5(g {0 \k))dk. 



(10) 



This solves the task of computing the spectrum of our example of the three-bond dressed 
star graph explicitly. 

In general, given a spectral equation (1) which does not fulfill (3), we generate a chain 
of derivative spectral equations g( m \k), where g( m \k) is the m'th derivative of (1) divided 
by S™, explicitly given by 



N 



g( m \k) = cos(S k — 7T7 + mix /2) — of 1 ' cos(Sjk — n^j + TO7r/2), 



(11) 



where af 1 ^ = af' (Sj/ So)" 1 . Since Sq < Sj, there always exists an M such that the amplitudes 



,(o) 



(m) 



satisfy the regularity condition (3), i.e., J2f=i < 1- Therefore, according to [7-9], 



root separators k^ exist on the level M and the roots k^ of g( M \k) = are explicitly 
computable via 

dg^ M \k) 



k (M) 



(M) 



k 



dk 



8{g {M \k))dk. 



(12) 



Since we now know the roots on the level M, we can go one step backwards to level M — 1. 
According to a root-counting argument [20] based on the Weyl formula [18,19], the root 
separators k^ 1 ' 1 ^ on the level M — 1 are the locations of the local extrema of g < ^ M ~ 1 \k) : 
which are given explicitly by the roots k^\ which we know. Therefore, fc^ M_1 ^ = k^ and 
the roots of g^ M ~ V) (k) = can now be computed explicitly, according to 



U(M-I) 



(M) 



kS n 



gM(k) 5{g^\k))dk. 



(13) 



Steps (12) and (13) define a recursive procedure, 

^ m " 1)= /(l, kS \g^{k)\S(g^- 1 \k))dk 1 m = M, M - 1, . . . , 2, (14) 



which can be followed until the level is reached and the roots k£\ i.e. the spectrum of the 
quantum graph, is known explicitly. 

It is important to notice that (12) - (14) are not just formal solutions. They yield 
m = 0, . . . ,M explicitly, by quadratures. Thus (12) - (14) constitute explicit solutions of 
the problem, very much in the spirit of the definition of explicit solutions by quadratures in 
the theory of differential equations [21]. 

Several special cases require discussion. If the regularity condition (3) is fulfilled, a root 
k n lies strictly inside of the interval [k n -i, k n ]. However, if (3) is not fulfilled, it is possible 
that a root k^ coincides with one of its separators k^} 1 or k^. This is, e.g., the case in 
our star-graph example above, where kf^ = n is a root and a root separator of (7). In the 
parameter space of a's and /3's cases like this are extremely rare (nongeneric). But even if 
such a case occurs, it does not present a problem for our theory. At the contrary, it saves 
one integration step since it can always be checked before performing the integration in (14), 
whether one of the separators k!^} 1 or is a root of g^ m ~ 1 \k). If so, the result k^L^ = 
k ( ^ 1 , or k^-^ = /4 m \ respectively, is known in advance, without actually performing the 
integration. 

In other special cases the roots of g( m \k) = can be obtained in the form of explicit 
periodic orbit expansions [7,8]. In order to illustrate this, let us return to our example of the 
three-bond star graph. We notice that the spectral equation g^(k) = of the three-bond 
star graph looks the same as the spectral equation [7] 

#4V-chain(*0 = sin^ofc) + ^2 sin( 1 S 1 /c) + r 2 r 3 sin{S 2 k) - r 3 sin{S 3 k) = (15) 

of the dressed four- vertex chain graph shown in Fig. lb, where r 2 = (f3\ — fa) /(Pi + P2), 
r 3 — (P2 — Pz)/(Pi + ^3) are the reflection coefficients at the vertices number 2 and 3 of 
the chain graph, and the actions So, . . . , S3 are the same as in (4). If we arrange for the 
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bond actions of the chain graph to equal the bond actions of the three-bond star graph, 
and furthermore arrange for = —r 2 , = —r 2 r 3 , = r 3 , which is possible if the 
scaling constants of the three-bond star graph fulfill (3f — fi\ + /?J = 0, then g^ l \k) of the 
three-bond star graph is the same as the spectral equation (15) of the associated four- vertex 
chain graph and the spectral points can be stated immediately and explicitly in the 
form of convergent, periodic orbit expansions [7,8], bypassing any integrations that would 
have been necessary according to the scheme defined in (14). 

Although they presented the first examples of explicitly solvable quantum graphs, a 
major shortcoming of [7] and [8] is the fact that the theory presented in [7] and [8] is 
only applicable to regular quantum graphs, i.e. quantum graphs that fulfill the regularity 
condition (3). In this paper we showed that the restriction to regular quantum graphs is 
not necessary: all scaling quantum graphs can be solved explicitly. Nevertheless, the theory 
presented in [7-9] provides an indispensable foundation without which the present theory 
would not be possible. 

A conceptual advance is the following. Frequently an operational definition of quantum 
chaos, or a quantum chaotic regime, is the "loss of quantum numbers" . To illustrate, let us 
consider a Hamiltonian system with Hamiltonian H = H + fiV, where H is an integrable 
Hamiltonian, \i is a real parameter and V, with respect to and in conjunction with H , is a 
nonintegrable perturbation. Many quantum systems, for instance the hydrogen atom in a 
strong magnetic field [10], can be described in this way. For /i = the system is integrable 
and possesses a complete set of quantum numbers 

that can be obtained, at least approximately, using EBK quantization [18,19]. As the 
parameter \i increases, EBK quantization breaks down and the system makes a transition to 
quantum chaos. This explains the frequently employed practice of characterizing the onset of 
quantum chaos by a loss of quantum numbers, since the breakdown of the EBK quantization 
scheme implies the loss of quantum numbers. The results obtained in this paper, however, 
show that this is not necessarily a good way to characterize quantum chaos. Although not 
strictly chaotic in the classical limit (due to ray-splitting [11,12,22,23] the term stochastic 
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may characterize the situation better), quantum graphs were shown by many authors [1-5] 
to be excellent models of quantum chaos. Yet, our results above show that a well-defined 
quantum number, the counting index n, still exists, and produces explicit energy levels in 
exactly the same spirit as the EBK quantization scheme. 

The iteration scheme (14) is perhaps the most interesting feature of our method of ex- 
plicitly solving quantum graphs. We call the smallest M that "regularizes" a given quantum 
graph (i.e. the amplitudes of g( M ^ fulfill (3)), the order of the quantum graph. For any given 
quantum graph its order is unique. Since the order M of a quantum graph determines the 
length of the bootstrapping iteration scheme (14), it is possible that the order of a quantum 
graph is also an indication of the complexity of its spectrum. Quantum iterations similar to 
(14) were studied before [24] and were found to lead to sensitivity and chaos on the quan- 
tum level. This may explain the reason why certain quantum graphs are such good models 
of quantum chaos [1-5] and the order M of the quantum graph may be an indication of 
how well a given quantum graph can be described in terms of the usual diagnostic tools of 
quantum chaos, such as, e.g., random matrix theory [18,19,25]. 

The authors acknowledge financial support by the National Science Foundation under 
Grant No. 9984075. 



8 



REFERENCES 

[1] T. Kottos and U. Smilansky, Phys. Rev. Lett. 79, 4794 (1997); Ann. Phys. (N.Y.) 274, 
76 (1999). 

[2] H. Schanz and U. Smilansky, Phys. Rev. Lett. 84, 1427 (2000). 
[3] T. Kottos and H. Schanz, Physica E 9, 523 (2001). 

[4] U. Smilansky, in Mesoscopic Quantum Physics, Les Houches, Session LXI, 1994, edited 
by E. Akkermans, G. Montambaux, J.-L. Pichard and J. Zinn- Justin (Elsevier Science 
Publishers, Amsterdam, 1995), pp. 373-433. 

[5] P. Kuchment, Waves in Rand. Media 12(4), Rl (2002). 

[6] Y. Dabaghian, R. V. Jensen, and R. Blumel, Phys. Rev. E 63, 066201 (2001). 

[7] Yu. Dabaghian, R. V. Jensen, and R. Blumel, JETP Lett. 74, 235 (2001); JETP 94, 
1201 (2002). 

[8] R. Blumel, Yu. Dabaghian, and R. V. Jensen, Phys. Rev. Lett. 88, 044101 (2002); Phys. 
Rev. E 65, 046222 (2002). 

[9] R. Blumel, Yu. Dabaghian, and R. V. Jensen, Mathematical foundations of regular 
quantum graphs, to be published (2003). 

[10] H. Friedrich, Theoretical Atomic Physics, second edition (Springer, Berlin, 1998). 

[11] R. E. Prange, E. Ott, T. M. Antonsen, B. Georgeot, and R. Blumel, Phys. Rev. E 53, 
207 (1996). 

[12] L. Sirko, P.M. Koch, and R. Blumel, Phys. Rev. Lett. 78, 2940 (1997). 
[13] M. Keeler and T. J. Morgan, Phys. Rev. Lett. 80, 5726 (1998). 

[14] Sz. Bauch, A. Bl§dowski, L. Sirko, P.M. Koch, and R. Blumel, Phys. Rev. E 57, 304 
(1998). 

9 



[15] R. Blumel, P. M. Koch, and L. Sirko, Found. Phys. 31, 269 (2001). 

[16] Y. C. Lai, C. Grebogi, R. Blumel, and M. Ding, Phys. Rev. A 45, 8284 (1992). 

[17] A. Kohler and R. Blumel, Phys. Lett. A 238, 271 (1998). 

[18] M. Gutzwiller, Chaos in Classical and Quantum Mechanics (Springer, New York, 1990). 

[19] H.-J. Stockmann, Quantum Chaos (Cambridge University Press, Cambridge, 1999). 

[20] Yu. Dabaghian and R. Blumel, in preparation for submittal to Phys. Rev. E (2003). 

[21] W. W. Stepanow, Lehrbuch der Differentialgleichungen (VEB Verlag, Berlin, 1976). 

[22] L. Couchman, E. Ott, and T. M. Antonsen, Jr., 
Phys. Rev. A 46, 6193 (1992). 

[23] R. Blumel, T. M. Antonsen, Jr., B. Georgeot, E. Ott, and R. E. Prange, Phys. Rev. 
Lett. 76, 2476 (1996); Phys. Rev. E 53, 3284 (1996). 

[24] R. Blumel, Phys. Rev. Lett. 73, 428 (1994). 

[25] F. Haake, Quantum Signatures of Chaos (Springer, Berlin, 1991). 



10 



Figure Captions 

Fig. 1: (a) Dressed three-bond star graph and (b) dressed four- vertex chain graph. Different 
potential strengths on the bonds are indicated by different thickness of the bonds. Different 
vertex strengths are indicated by different dot sizes representing the vertices. 



11 



